Can Cosmic Rays Heat the Intergalactic Medium? 
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Supernova explosions in the early star forming galaxies will accelerate cosmic rays (CRs) . CRs are 
typically confined in the collapsed objects for a short period before escaping into the intergalactic 
medium (IGM). Galactic outflows can facilitate this escape by advecting CRs into the IGM. An 
outflow that results in a termination shock can also generate more CRs. We show that the CR 
protons from the above processes can significantly affect the thermal history of the IGM. Within 
plausible range of parameters, cosmic ray heating can compensate for adiabatic cooling and explain 
the measured IGM temperature at redshifts z ~ 2 — 4, even with early reionization. 

PACS numbers: 98.62.Ra, 98.70.Sa, 96.40-z, 98.62.-g, 98.80.-k 



It is widely believed that the temperature of the low 
density intergalactic medium retains some memory of the 
reionization process [jj] . Analysis of Lya absorption spec- 
tra in the framework of hydrodynamical simulations and 
semi- analytic models |2j suggest that, the mean temper- 
ature of the IGM is in the range 10 4 < T(K) < 3 x 10 4 
for 2 < z < 4. This has been used to argue that hy- 
drogen reionization occurred below z ~ 9 3]. Had the 
reionization occurred at higher z, adiabatic cooling would 
have brought the temperature much below the measured 
value at z = 4. The detection of a strong Gunn-Peterson 
trough and the sizes of ionized regions around the highest 
z QSOs 4], arc consistent with a large neutral hydrogen 
fraction at z ~ 6. The CMB data from the WMAP satel- 
lite 0], however, points to an earlier epoch of reioniza- 
tion in the range 11 < z < 30. Thus either new heating 
sources (other than photo-heating) or multiple episodes 
of reionization with heating at low z provided by, say, 
the He II reionization induced by QSOs 0, are needed 
to explain the IGM temperature at z ~ 2 — 4. Here, we 
consider the effect of cosmic ray protons on the thermal 
history of the IGM. 

Cosmic rays are expected to be generated and acceler- 
ated from the shocks created by the exploding supernovae 
(SNe) in the star forming regions, termination shocks cre- 
ated by outflows from galaxies and accretion shocks dur- 
ing structure formation. The properties of CRs in our 
Galaxy is well documented 0, 0, E3 • The energy den- 
sity in the proton component is about 1 eV cm -3 , and 
CRs are thought to be confined to the Galactic disk for 
about 10 7 yr before escaping, presumably into the IGM. 
It is widely accepted that > 15% of the average SNe en- 
ergy must go into accelerating the CRs so that the flux 
density of CRs can be maintained at the observed value 
|lflj . Low energy CR protons play an important role in 
the ionization and thermal state of the gas in the inter- 
stellar medium (ISM) of normal galaxies. The possibility 
that they could also play an important role in determin- 
ing the thermal history of the IGM has been considered 
sporadically in the literature 0, ■ The existence of 
good observational constraints on (i) the global star for- 



mation rate (SFR) density as a function of z (ii) epoch 
of reionization and (iii) temperature of the IGM makes it 
appropriate to re-examine this issue using a model that 
self consistently computes the generation and spectral 
evolution of CRs in the IGM. We do this here in the 
framework of a standard flat ACDM model (fi m = 0.27, 
n b = 0.04 and H = 70 km s _1 Mpc" 1 ). 

Consider CR protons produced from the SNe explo- 
sions. To compute the rate of energy injected into CR 
protons from early star forming galaxies, one needs to 
know the SFR and the number of SNe explosions, fsN, 
per solar mass of forming stars. We use the observed 
global comoving SFR density [13, i?(z) measured in units 
of Mq yr -1 Mpc~ 3 and approximated by 
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r), in the range 5 to 9, is the unknown sc aling 
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factor used to correct the effect of dust reddening |14| . 
We shall adopt rj = 7 in our models. Further, using 
the "Starburst99" code GJ] we get f SN = 8x 10~ 3 for 
Salpeter IMF (with M min = 0.1 M , M max = 100 M ) 
and fsN — 1.7 x 10~ 2 for a top-hat mass function with 
M min = 40 M , M max = 100 M (denoted as "top- 
heavy"). We assume that the energy released from a 
single SNe is 10 51 erg and a fraction e = 0.15 of this is 
utilized for accelerating the CR protons. The average 
rate of energy injected in the CR protons is then, 



E CR {z) = 10- 30 e R f SN (l + zf erg s" 



1 cm 3 . 



(2) 



The cosmic ray proton spectrum from the SNe is in gen- 
eral taken to be a power law in momentum space, with 
n(p) dp oc p~ q dp, with q — 2.2 as a typical value Q. We 
can then compute the average rate of CRs generated in 
star forming regions as 



dn(p, z) 
Jt 



dp = N p~ q dp 



(3) 



where the normalization Nq for a given z is fixed by inte- 
grating E(p)N p~ q dp using a lower energy cut-off Q of 
10 keV, and equating the result to Ecr(z). Here E(p) 
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FIG. 1: Cosmic rays spectrum originating from SNe (dot- 
ted one) and the modified spectra (solid lines) escaping from 
galaxy. The curves A and B are for x — 10 and 3 g cm -2 
respectively. The curve C is for the advection model, where 
the spectrum A has been modified by decreasing the momen- 
tum of each CR particle by a factor 10. Curve D is the CR 
spectrum originating from the termination shock with a sharp 
cut-off at 3 MeV. For illustrative purpose the injected Q(E, z) 
is normalized using R at z — 15. 



is the kinetic energy of a CR proton corresponding to a 
momentum p. This injected spectrum will be modified 
as the CRs propagate first within the high redshift star 
forming galaxy and then into the IGM. The evolution 
of the velocity v = j3c of the CR particle due to colli- 
sions with free electrons, hydrogen ionization, and the 
adiabatic expansion of the universe is given by [21 Il2| , 



^ = 3.27 x 10- 16 n e (z) 
dt K ' 



(1-/32)3/2 p 



1.94 x lCr 16 n ff/ (z) 
2f3 2 



(3 a& + Z? 3 

(1-/32)3/2 
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[(0.01) 3 + 2/^ 
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[1 + 0.0185 ln/3 6(/3-0.01)] 
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dt 



Here, n e and njjj are the densities of electrons and neu- 
tral hydrogen respectively in units of cm~ 3 . is Heavi- 
side function, x rn = 0.0286 (T/2 x 10 6 K) 1 / 2 and T is the 
temperature of the gas. The energy losses are dominated 
by low energy CRs for 3 > x m . Typically a CR parti- 
cle of energy E will lose all its energy to the electrons in 
a time, t L ~ 1.5 x 10 8 (E/3MeV) 1 - 5 (W- 4 '/n e ) yr, after 
traveling a distance d L ~ 4(-E/3MeV) 2 (10~ 4 / ri e) Mpc. 

The escape of the CRs from the star forming galaxies 
will depend on the structure and strength of its magnetic 
field, the gas density (p) and outflows etc. We consider 
three possible scenarios for the emerging CR spectrum: 
(I) CRs are confined for a period t c to the galaxy and then 
escape into the IGM, (II) after a confinement time t c , 



CRs are advected by galactic outflows into the IGM, and 
(III) CRs are generated in termination shocks created by 
the outflows. 

In case (I), we compute the emergent spectrum from 
the galaxy, taking into account the energy loss to elec- 
trons, as a function of grammage, x ~ pct c g cm~ 2 , tra- 
versed by the CRs within the galaxy. For our Galaxy 
x is ~ 10 g cm -2 |9j. For the high z galaxies, x could 
be smaller, because the global star formation rate den- 
sity is expected to be dominated by smaller mass ob- 
jects that may also be 'magnetically' younger and will 
confine CRs for a smaller period. The injected spec- 
trum as well as modified spectra (defined as the rate 
Q(E,z)dE = (dn/dt)(dt/dz)dp per unit redshift) for 
x = 10, 3 g cm -2 , are shown in the Fig. H an d labeled as 
A and B. Larger and larger part of the low energy CR- 
spectrum is cut off due to energy losses as one increases 
x. As an example, when x — 10 g cm" 2 protons of initial 
energy < 160 MeV will lose all their energy. Energy loss 
from protons with higher energy does lead to an universal 
low energy tail with Q oc E x l 2 . Since dE/dt oc i?~ 1//2 , 
the heating is still dominated by protons near the peak 
of Q. 

High redshift Lyman bre ak g alaxies (LBGs) show sig- 
natures of strong outflows [l6j. Such outflows, that are 
also required to explain the observed IGM metallicity, 
will help in advecting the CRs into the IGM. The adia- 
batic loss in such outflows can also lower the momentum 
(p) of the CR particles, with p oc p 1 ^ 3 , thereby shifting 
the spectrum to lower momentum and decreasing the to- 
tal energy in CRs. Suppose a wind from a galaxy of 3 
kpc size expands up to 30 kpc before CRs escape into the 
IGM. The average gas density will decrease by a factor 
10 3 there by reducing the momentum of CRs by a factor 
10. The emergent spectrum for this case (our case (II)) 
(and x = 10 g cm -2 ) is also shown in Fig. (see curve 
C). 

In order for the CRs to influence the properties of the 
IGM, they should be released into the IGM within the 
Hubble time. This time-scale is confinement time for case 
(I) and advection time for case (II). In our Galaxy t c ~ 
10 7 yr in the disk. In the Galactic halo of say a radius 
Th, the typical t c is the diffusion time, r 2 H /QD ~ 10 s yr, 
for th ~ 10 kpc and a diffusion coefficient D = 5 x 10 28 
cm 2 s _1 |8J. Thus for case (I) t c is expected to be smaller 
than the Hubble time, even for the Galaxy. Similarly the 
advection time scale in an outflow of velocity typically 
seen in LBGs (i.e ^100 km s _1 ) to a radius 30 kpc as in 
our case (II) will be 3 x 10 8 yr, which is shorter than the 
Hubble time for z < 15. 

In case (III) we consider the CRs produced from the 
termination shocks of the outflows. We assume that 
about 10% of the SNe energy is channeled into driving 
such an outflow, and a fraction e of this goes into accel- 
erating the CRs in the resulting termination shock. The 
spectrum of CRs emerging from the shock is again ex- 
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pected to be a power law in momentum space. As the 
wind termination shock is expected to form far out in the 
halo of the galaxy and propagate for most periods in the 
IGM itself, even low energy cosmic rays can be directly 
injected into the IGM. However the lowest energy CRs 
will lose their energy in the vicinity of the galaxy itself. 
For example, a 3 MeV CR proton will travel a proper 
distance g?l ~ 10 Mpc in the IGM at z ~ 5 before losing 
all its energy. While within a halo that has 200 times 
larger density, it will travel ~ 50 kpc. Thus we assume 
that the particles with energy E > 3 MeV will be able 
to escape from termination shocks and contribute to the 
heating of the IGM. The resulting spectrum with a sharp 
cut-off at 3 MeV is shown in Fig. \T\as case D. 

The free propagation of CRs in the IGM will depend on 
the strength and coherence of the intergalactic magnetic 
field. For example in the case of ISM the streaming of 
CRs is limited by resonant scattering on self-generated 
Alfvcn waves. It turns out that a coherent IGM magnetic 
field larger than ~ 10~ 14 G is required for a streaming 
instability growth time-scale < 10 8 yr (Eq. 32 in 0|), so 
that there can be adequate growth of Alfven waves, by 
say z ~ 5. Astrophysical battery mechanisms produce 
only tiny coherent IGM seed fields ~ 10~ 23 - 10~ 19 G 
17]. Early Universe generation mechanisms are not very 
compelling because the results are often exponentially 
sensitive to parameters involved ^^]. Further, for z > 2, 
the self-confinement due to plasma waves generated by 
the streaming CRs ^l) is inefficient because of the strong 
damping of waves by electron ion collisions. Here, we will 
assume that the CRs escaping the high z galaxies are not 
self confined, but can stream out freely in the IGM. 

The spectrum of cosmic ray background at a given z 
is the sum of CR spectrum injected at z and that evolv- 
ing from higher redshifts. The number density of CRs, 
Ncr{/3, z, zo)df3, in an interval (/3, j3 + d/3), at any given 
epoch, z, is 



Ncr(P, z,z ) 



dz 



dn(zj,pi) dpi d(3j ( 1 + z 
dzi dj3i df3 \ 1 + z% 



(5) 

In order to compute Nqr one needs to know the /? evo- 
lution and hence the ionization history of the IGM. We 
model this using a very simple rcionization scenario with 
the aim to merely compute Ncr in a self-consistent set- 
ting rather than addressing the issues related to early 
rcionization. The rate of change in the number density 
of ionized hydrogen (i.e. nn n) can be written as, 



dt 



A 7 + / N CR ((3, z, z )cn H ia H ij3d(3 



a B n HII n e . 



(6) 



The right hand side of Eq. H3 gives the rates of photoion- 
ization, cosmic ray ionization and recombination respec- 
tively. We assume all the Lyman continuum photons 




FIG. 2: Predicted temperature evolution of the IGM for dif- 
ferent injected spectra given in Fig. We follow the labeling 
convention as in Fig. Temperature evolution for the adia- 
batic cooling is shown with a dotted curve. The data points 
are a compendium of observations presented in 0. 



entering the IGM with a rate N 7 are used for ionizing 
H I. For a SFR of 1 M yr" 1 about 24000 s" 1 and 4000 
s _1 Lyman continuum photons are generated per baryon, 
for "top-heavy" and "Salpeter" IMF respectively. This 
is used with an escape fraction of 0.1 to compute an 
is the cross-section for H I ionization by protons [20J and 
as is the case-B recombination coefficient. We find CR 
ionization is sub-dominant compared to the photoioniza- 
tion. 

We assume "top-heavy" IMF prior to the epoch of 
rcionization and "Salpeter" IMF after that. The value 
of ~a.Hu at a given z is obtained using Eq. [S] and the an- 
alytical fit to the observed SFR density as a function of 
z (Eq. ^ . We get the epoch of reionization (where uhii 
— tih) at z r ~ 10. Note that this is independent of our 
choice of r\ because while R cx r\ the UV escape fraction 
f e cx I/77. However, as SNe rate increases with n the 
CR flux density will depend on rj. Using the reionization 
model, and Eqs. QJ-^wc compute the spectrum of CR at 
any given z, starting at zq = 15. Such spectra are used 
below to compute the heating of the IGM due to CRs. 

In our model calculations we assume (i) IGM reaches 



an uniform temperature of T = 3 x 10 K at z 



(A) 



heating by UV photons, Compton cooling by CMB and 
recombination cooling are negligible for z < z r and (iii) 
for simplicity the IGM is assumed to be of uniform den- 
sity. Thus the thermal state of the IGM after reionization 
is described by, 



d£ 

dz 



58 



F 



CR 



(7) 



with £ is the thermal energy per unit volume of the IGM. 
The first term in the right hand side of Eq. gives adi- 
abatic cooling rate. The second term Tcr. is the cosmic 
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rays heating rate given by 



dz 

dE(/3) dt 



dt dz 



N CR (f3,z,z )dp 
N CR {P,z,z )d(3. 



(8) 



Here, dE/dt is obtained from the first term of Eq. 0] 
Using £ — 3nkT/2, where T is the temperature and n the 
number density of particles in the IGM, one can obtain 



strained by observations and plausible range of parame- 
ters, the observed temperature of the IGM at 2 < z < 4 
can be explained as due to the CR heating, even if the 
epoch of rcionization is as early as z = 10. This implies 
that low energy CRs in the IGM can erase the thermal 
memory of reionization. In conclusion, cosmic rays do 
provide an important source for heating the intergalactic 
medium, that can explain the measured IGM tempera- 
tures at high redshifts. 



T(z) = T(z r ) 



1 + z 
1 + z r 

(l+Z r ) 2 

3n kT(z r ) J (1 + z'f 



CR 



(z',z ) 



dz' 
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(9) 



where tlq = n e (z = 0). 

The temperature evolution of the IGM for different 
injected spectra (shown in Fig. ^| computed using Eq. El 
are shown in Fig. |2 The temperature of the IGM goes 
well below the observed range due to adiabatic cooling in 
the absence of heating sources (dotted curve in Fig. |5J). 
However for all the models discussed here, CR heating 
leads to, T > 10 4 K in the rcdshift range 2 < z < 4. 

For case (I) the average grammage traversed by the 
CRs in the high z star forming galaxies has to be less 
than 10 g cm -2 seen in the Galaxy, in order to produce 
the observed mean IGM temperature (curves A and B 
in Fig. [21 • A lower grammage is reasonable for the high 
z galaxies, as mentioned earlier. These models predict 
cosmic ray energy density U cr ~ 6 x 10~ 3 eV cm~ 3 at z = 
0. Models with adiabatic losses (case (II)) or termination 
shocks (case (III)) produce higher temperatures, as they 
efficiently inject lower energy CRs into the IGM. At the 
same time they lead to overall lower cosmic ray energy 
densities of the order 1.5 X 10~ 4 eV cm -3 and 5.7 x 10 -4 
eV cm~ 3 for case (II) and (III) respectively at z — 0. 

Models of shock acceleration [lOj predict e as high as 
0.3 — 0.4. In such a case even a model with grammage 
as seen in our galaxy will produce higher temperatures 
than the observed mean. Simulations of pregalactic out- 
flows [2l]], find up to 30% of the supernovae energy can 
go into winds. The resulting temperature in case (III) 
can thus be much higher. We considered CRs produced 
in the SNe and termination shocks, CRs generated in 
structure formation shocks |22j could provide additional 
heating. The major uncertainties are related to the ex- 
tent to which the IGM is magnetized and confinement 
issues. More detailed computations in the framework of 
structure formation models are in progress. 

We have explored the influence of CRs escaping from 
the early star forming regions on the thermal history of 
the IGM. The estimated IGM temperature evolution is 
sensitive to the star formation rate (rf) and the confine- 
ment of CRs. Using a star formation rate density con- 
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